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I. INTRODUCTION 



The detection of gravitational waves is arguable one the most important outstanding exper- 
imental challenges in physics. With the construction of laser interferometric gravitational wave 
detectors like LIGO, VIRGO, TAMA, GEO600 there are good chances of direct detection in the 
very near future fl, Q, fl, fl, fl, fl, fl. Along side the interferometers, which are mainly 
aimed at detecting gravitational waves of astrophysical origin, the anisotropies in temperature and 
polarization of the Cosmic Microwave Background (CMB) have a strong potential to discover relic 



gravitational waves 0,0], [3] (see 



12l ] for a recent reviews). 
Most of the current techniques to detect gravitational waves are based on their interaction 
with electromagnetic fields. In general, the interaction of gravitational waves with electromagnetic 
radiation leaves imprints on the latter that can be experimentally measured [l^. In this work we 
shall deal with one such interaction effect which we shall call the "surfing effect" , where (figuratively 
speaking) the electromagnetic wave surfs on a gravitational wave leading to an observable phase 



change in the electromagnetic wave. This effect was first considered in 14|. In the present paper, we 
shall expand on the results of [l^ , [3] and generalize the effect for an arbitrary gravitational wave 
field. We shall consider the consequences of this effect for the planned precision radio (or x-ray) 
interferometric projects [3], [l^]. As we shall show, due to the transverse nature of gravitational 
waves, the surfing effect leads to an observable phase change only when the velocity of gravitational 
waves is different from speed of light. Using this fact, and the micro- arcsecond accuracy promised 



171 ]. we can place significant upper limits on 



by the precision interferometry measurements [l^ . 
the parameter e = (c — fg«,)/c which characterizes the deviation of velocity of gravitational waves 
from speed of light. 

The constraints on the speed of gravitational waves is an interesting experimental challenge. 
A potentially strong method of constraining e is to compare the arrival times of a gravitational 
wave and an electromagnetic wave emitted by a supernova or a gamma ray burst [l8(] (see also 



19[). Although this method will be able to give very strong constraints, it crucially depends on our 
ability to model and detect the gravitational wave signal from these sources, which is a significant 
theoretical and experimental challenge. In [3] the speed of gravity, and correspondingly indirectly 
the speed of gravitational waves, was constrained, by analysis of retarded gravitational potentials 
in the non-wave zone, by measuring propagation of the quasar's radio signal past Jupiter. Let us 
note, without going into detail of ambiguity in the interpretation of this result (see for example 



2ll]). that speed of gravity was constrained to Cg = (1.06 ± 0.21)c corresponding to e < 0.2. In this 
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paper we consider the propagation of electromagnetic radiation in the field of gravitational waves, 
i.e. in the wave zone of the gravitational field, where the velocity parameter e can be introduced 
avoiding any ambiguity. Using the surfing effect, we consider an independent method of placing 
upper limits on the e-parameter, which could potentially give very strong limits on the speed of 
gravitational waves (e < 10^^ for 10^-'^''). 

The plan of the paper is as follows. In Section [II] we shall consider the propagation of an 
electromagnetic wave in the field of a single monochromatic gravitational wave. We shall discuss 
the physical aspects of the surfing effect with a view on the precision interferometric measurements, 
as well as write down some of the equations that will be used in the following sections. In Section 
mil we generalize the surfing effect for an arbitrary gravitational wave field. Positing statistical 
properties of the gravitational wave field we derive the consequent statistical properties of the 
response of an interferometer. In Section IIVI we use the surfing effect along with the predicted 
precision level of the interferometry measurements to place upper limits on the velocity parameter 
e depending on the energy density of gravitational wave described by the density parameter 
Finally, in Section |V] we present a short discussion and summary of the main results of this paper. 



II. SINGLE MONOCHROMATIC GRAVITATIONAL WAVE 



Let us consider a slightly perturbed flat Priedmann-Lemaitre-Robertson- Walker (FLRW) uni- 



verse with coordinates = {ri,x^) and the metric given by 
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2M 



ds^ = g^^udx^dx" = o?{-r]) [-d-rf + {5ij + hij) dx^dx^] . (1) 

Here rj is conformal time coordinate related to proper time coordinate t through a relation cdt = 
a{rj)dr], and are the spatial coordinates. In the above expression a(r/) is the scale factor, and hij 
is the gravitational wave perturbation. Since, in this work we shall be interested in electromagnetic 
waves emitted by objects with a redshift z < 5, we will work with a matter dominated Universe 
model characterized by a scale factor 

a(r?)=a„(^y, (2) 

where rjo corresponds to present time. From the above expression it follows that is the scale 
factor at the present time. Without loss of generality we can set Oq = 1. With this convention, the 
value of ijo is related with the present day Hubble constant Hq 

Vo = ^ = '2Lh, (3) 



where we have introduced the Hubble length Lh = c/Hq. 

For simplicity, in this section we shall consider the case of a single monochromatic plane grav- 
itational wave. We shall also restrict our considerations to gravitational waves whose wavelength 
at the current epoch is small compared to the present Hubble length. This implies that we can 
consider these waves to be monochromatic waves with an amplitude damping adiabatically with 
the expansion of the universe 

where is the amplitude at the present epoch, is the polarization tensor of the gravitational 
wave, and = (ko, hi) is the wave vector of the gravitational wave 



22 . It is convenient to 
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introduce the wavenumber k = [dijk'^k^^ . The present day wavelength of the gravitational field 
Xgw, is related to wavenumber k by the relation A^^ = 27: /k. The present day frequency of the 
gravitational wave tOgw is related to the time component of the wavevector /cq by relation a;^^ = cko . 

Let us analyze the electromagnetic wave propagation in the approximation of geometrical optics. 
In this approximation, the wave equation is written for the quantity ip{x^) = ^{r],x^), known as 



the eikonal [24]. Eikonal has the physical meaning of the phase of an electromagnetic wave field 
f(x^) = A{x^)e'''^^^'^\ (The quantity A(x^) describes the amplitude of the wave field, but shall 
not interest us in this consideration.) The wave vector of the electromagnetic wave is give by 

The eikonal equation follows from the isotropy condition = 0. Substituting the expres- 



sion for in terms of if) we arrive at the eikonal equation 
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We shall seek the solution to equation dH) in a perturbative form 

i,{x^) = i,Q{x^)-r^x{x^). (6) 

Here the zeroth order solution corresponds to the solution in the absence of perturbations, while 
the first order solution corresponds to the solution in the presence of perturbations ~ h. Let us 
assume that the zeroth order solution corresponds to a plane monochromatic electromagnetic wave: 

-(Ao = WE (f/ + eiX*) /c, 

(7) 

where e* {S^^eiCj = 1) is the unit vector along the direction of wave propagation, and a;^; is the 
frequency of the electromagnetic wave at the present time. Taking into account the solution ([7]) of 



the zeroth order equation, the first order equation for ipi^x'^) takes the form 

The solution to equation ^ is written in terms of the hne of sight integral (along the unperturbed 
light path): 

2c 



il^ii'n.x') = ^ujE e'e^ I ds hik{s,x' + e' {s -rjo)) , (9) 



Vo-D 

where s is the j^-time parameter along the light ray path from the emitter to observer. The limits 
of integration (rjo — D) and rjo correspond to the time of emission and observation correspondingly. 
Using ([2j), the parameter D can be related to the redshift of emission z by 

" = - TiTi) ■ <^°' 

Below, for simplicity of analysis and in order not to obscure the physical interpretation of 
the surfing effect, we shall consider the problem in flat space-time without cosmological evolution 
of the scale factor, and correspondingly, assume no cosmological evolution of gravitational wave 
amplitude. This analysis is equivalent to setting the scale factor in expression ([1]) to a constant 
value, i.e. a['q) = = 1. Alternatively, this approximation can be viewed as an analysis restricted 
to small values of redshift, i.e. z <^1. In this limit, parameter D represents the physical distance to 
the source, and is related to the redshift by the usual Hubble law D ~ zLh- We shall reintroduce 
the cosmological evolution in Section IIIIl where we shall explain how and why the result modifies. 
Some of the calculational subtleties that arise when analyzing this situation are considered in 
Appendix [B1 

Thus, assuming a plane gravitational wave with constant amplitude (i.e. setting a{r]) = ao 
in ([4])), we can explicitly evaluate the integral in ([9]) to get the resulting phase change due to a 
gravitational wave 

7 / JD(ko-he^) i\ 

Mv,x^) = ^ e^eV.. e-C^^^-'^-O ^ , , • (11) 

2c y ^{Ko - he") j 

An interferometer is an experimental device capable of measuring the difference in phase in an 
electromagnetic wave jlT^ . The micro-arcsecond precision promised by the upcoming interferometry 
projects may allow to detect signature of gravitational waves or set upper limits on their magnitude. 
For this reason, let us switch our attention to the interferometers. The interferometers measure the 
variation in the phase of the electromagnetic signal received by its base antennae. Let us consider 



a long base interferometer setup. In this case there are two antennae separated by a spatial vector 
U. Without loosing generality, we can assume that the first of this antennae is located at the 
coordinate origin. Then, the difference of phase measured by the by the two antennae, due to the 
gravitational wave influence, is given by 

AV^i = V'i(r/,0)-^i(7?,L^). (12) 



For the plane electromagnetic wave under consideration, substituting (llip into (|12p we arrive at 
the expression for phase difference measured by the interferometer 

gjD(fco-fcie') _ -j^ 



zc 



l-e*-^M . (13) 



The output of an interferometric measurements is usually quoted in terms of angular resolution. 
We shall use this convention. The angular resolution of an interferometer corresponding to a phase 
resolution l^ipi is give by 

Aa = -^AV'i, (14) 
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where L = (5ijUU^ is the baseline length of the interferometer. 

In what follows, we shall be interested in gravitational waves with wavelengths considerably 
larger that the baseline of the interferometer, i.e. L/Xg^ <C 1. In the case of a space-borne 
interferometer with base length of 10® /cm, this implies Ugw ^ 10^^ Hz for the gravitational wave 
frequency. For interferometers with a shorter baseline this limit can be significantly larger. 

To proceed, let us further introduce unit vectors = U /L and A;* = /k. is a unit vector 
pointing in the direction of the interferometer baseline, while A:* is the unit vector pointing in the 
direction of the gravitational wave propagation. Assuming L/Xg^ ^ 1 in equation (jl3p . we can 
expand this expression in a series retaining only the lowest order term in kiU. We get 

{ko - kiC') 

Up to now we have not posited any relationship between the gravitational wave frequency c/cq 
and the magnitude of the wavenumber k. This relationship (dispersion relation) defines the velocity 
of a gravitational wave Vgw = cko/k. In General Relativity this velocity equals the speed of light, 
but in an alternative theory this might not be the case. In order to analyze the possibility that 
Vgw 7^ c, let us use the phenomenological parameter e, introduced in the previous section 

C — Vguj , Cko ,^ „K 

e = where Vgy, = —, (16) 



Aq = ^ho krF e^e^Pik e"^^"^ 



(15) 



(18) 



which characterizes the relative deviation of velocity of gravitational waves from the speed of light. 
Let us note that, e can be related to a non vanishing rest mass of a graviton rUg through the 
relation 

e = l- ^ (17) 
Returning to equation p^ . and substituting the relationship = {1 — e)k into it, we get 

1 — e — fcjC* 

It is instructive at this point to look more closely at expression (jlSp . We are considering an elec- 
tromagnetic wave travelling "along" a plane gravitational wave. As follows from this expression, 
the angular displacement |Aa| is most pronounced when both the waves move in almost parallel 
directions. This picture is reminiscent of wave surfing, and hence we call this the "surfing" ef- 
fect. The expression in square brackets becomes large (proportional to kD ~ D/Xg^), when the 
denominator tends to zero (i.e. 1 — e — /cje* — > 0), leading to a resonance effect. In the case when 
e = 0, this does not lead to significant growth of the angular displacement |Aa|, because of the 
transverse nature of gravitational waves (since e^e^pij ^ as /cje* 1). On the other hand, if 
e 7^ 0, the expression for |Aa| becomes sufficiently large for a gravitational wave propagating at 
an angle cos 9 = kie^ ~ (1 — e) to the line of sight. Thus, due to the transverse nature of the 
gravitational waves this surfing effect is absent if e = 0, but can become significant for e 7^ 0. This 
effect, as we shall show in the following section, can be used to put stringent constraints on the e 
parameter characterizing the velocity of gravitational waves. 

III. ARBITRARY GRAVITATIONAL WAVE FIELD 

In the previous section we considered the case of a single monochromatic gravitational wave. 
In order to generalize the considerations of the previous section, let us now consider an arbitrary 
gravitational wave field. This field can be decomposed into spatial Fourier modes 

h,{r^y) = I S\^Y. [hs{k\r^) (fe')e*-^ + K{k\r^) (A;')e-*-^] , (19) 

.8=1,2 

where d^k denotes the integration over all possible wave vectors, and s = 1,2 correspond to the 
two linearly independent polarizations of a gravitational wave. The mode functions hs{k^,ri) have 
the following time evolution 

hsik\v) = hs{k') e-^^(^-^)^ (20) 



where hs{k'^) is the gravitational wave amplitude at the present time. 

Due to the linear nature of the problem, following the decomposition (|19p . the total angular 
displacement due to gravitational waves can also be decomposed into Fourier modes in a similar 
fashion 

Aa= f d\ ^ [hsik')Aasik') + K{k')Aa*{k')] , (21) 

s=l,2 

where we have introduced a tilde over a, in the right hand side of the above expression, to indicate 
the explicit factoring out of the gravitational wave amplitude hs{k^) compared with expression 

dlED. 

Using the results of previous subsection (see Eq. (jlSh ). ignoring cosmological evolution for now, 
the contribution from a single Fourier component Q!s(/c*) is given by 

1 _ ^ikD{l-e-kie'] 



1 — e — he^ 



(22) 



In general for an arbitrary gravitational wave field ()19p the angular displacement measured by 
the interferometer is given by expressions (j21|) and (j22|) . In practice we do not have the precise 
information about the gravitational wave field, and are restricted to knowledge of only its statistical 
properties. Let us assume the following statistical properties for the mode functions hs{k^) 

<hs{k')>=0, <hs{k') h*Ak")>=^^6ss'6^{k' -k''), (23) 

where the brackets denote ensemble averaging, and Ph{k) is the metric power spectrum per loga- 
rithmic interval of k. These conditions correspond to a stationary statistically homogeneous and 
isotropic gravitational wave field. 

The statistical properties of Aa follow from the statistical properties of the underlying gravita- 
tional wave field ([23l) . Using ([2T]) , ([22]) and (f23]) after straight forward calculations we get following 
statistical properties for angular displacement Aa: 

< Aa > = 0, (24a) 
—Phik)Aa\k), (24b) 

where we have introduced the transfer function 

Aa'^ik) = — [ dnS" \Aas{k')f . (25) 
SttJ Y 

In the above expression represents integration over the possible directions of g.w. wave 
(i.e. d^k = k'^dkdQ). 




FIG. 1: The graphical representation of the various vectors and angles used in the text. Vector e* is the 
unit vector along the direction of the electromagnetic wave, unit vector P is aligned with the base of the 
interferometer, and fc' is the wave vector of the gravitational wave. 

We shall now proceed to calculate the expression (j25p explicitly. Let us introduce a spherical 
coordinate system (9, (p) related to the spatial coordinates {x*} in the usual manner 2J]. Without 



loss of generality, we can assume that we are looking in the north-pole direction, i.e. = (0,0, 1). 
Let us also introduce the quantity fi = cos 9 = eik^, characterizing the angle between a gravitational 
wave and the direction of observation. Furthermore for an interferometer, for optimal resolution 
is aligned perpendicular to e*, thus without loss of generality we can assume P = (1,0,0). The 
geometry of the problem is presented in Figured! The polarization tensors for gravitational waves 
have the form Pij (/c*) = (ef it ief ){e^ it ie'j )/2, with it corresponding to the two independent 
circularly polarized degrees of freedom s = 1,2 {ef and ef are the meridian and azimuthal unit 
vectors perpendicular to the gravitational wave wavevector ki, for a detailed discussion see for 



example 



25(1). Taking into account the relations 



and substituting (j22p into (j25p . after straight forward manipulations, the expression (j25p for the 
transfer function takes the form: 



Aa'^ik) 



1 

16 



+1 



fi/i (l — 



2\3 



-1 



sm' 



:i 



(26) 



The terms in the above integral have a clear physical meaning. The factor (1 — fi'^)^ is due to the 
transverse nature of the gravitational waves and the geometry of space interferometry ((1 — fj?)"^ 
and (1 — /x^) terms correspondingly). The quantity in square brackets sharply peaks at values 



/i ~ (1 — e), which is the result of a resonance effect, i.e. what we call the surfing effect, for 
gravitational waves travelling at an angle cos9 ~ (1 — e) to the line of sight. Due to the pre-factor 
(1 — ij?)^, this resonance does not give a significant contribution for the case e = 0. The integrand 
in expression (j26p is plotted for the two cases in Figure [2j In the limit e ^ and D/\gw ^ oo we 
can calculate the transfer function (j26p explicitly. Referring the reader to appendix |X] for details 
of calculation, let us present the final result below: 



1 

20 



1 + h-KC'kD 



(27) 



As was mentioned previously, when deriving expression (j27p we had ignored the cosmological 
evolution of the gravitational wave amplitude. In Appendix [B] we derive the expression for the 
transfer function when the cosmological evolution of gravitational waves is properly taken into 
account. The resulting expression is as follows: 



1 

20 
1 

20 



+ b-Ke-'kD (1 + z) 
+ bTre^kLn (1 + z) ( 2 



(28) 



From expression (j28p we can quantify the condition for the resonance to occur by comparing the 
two terms in the square brackets. The resonance occurs when, in the right side of (j28p . the second 
term (resonance term) is larger that the first term (non-resonance term), i.e. when e^kLu 3> 1. 
As we shall show in the next section, given the planned level of sensitivity for interferometric 
measurements, the resonance effect allows to place significant upper bounds on the parameter e. 

The redshift factors z, occurring in expression (I28p . have a clear physical interpretation. The 
factor (1 + z)"^ in the non-resonance part of the transfer function occurs because, its main con- 
tribution comes from epoch when the electromagnetic radiation was emitted, corresponding to a 
redshift of z. It is thus sensitive to the gravitational wave power spectrum at the epoch of emission, 
which was a factor (1 + z)^ stronger than today. On the other hand, in the case of the resonance 
term, the contribution to the transfer function is gained along the path from emitter to observer. 
This leads to a weaker z dependence in this term compared with the non-resonance term. 



IV. UPPER LIMIT ON THE VELOCITY OF GRAVITATIONAL WAVES 

Let us now consider the implications of the surfing effect for the precision interferometry mea- 
surements, and the achievable upper limits on e. When considering stochastic gravitational wave 
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FIG. 2: The illustration of the resonance effect, present for e 7^ 0. The graphs show integrand in expression 
(PS)) . For the case e 7^ the integrand sharply peaks at angle /i « (1 — e) (solid red line), while for the case 
e = the effect is absent (dashed blue line). In the case of e 7^ 0, the gravitational waves travelling at an 
angle cos 9 =« (1 — e) to the line of sight are the predominant contributors to the surfing effect. The figure 
on the left shows the integrand for the whole region of /i, while the figure on the right zooms into the region 
around the resonance. 

fields, it is customary to introduce the density parameter 0,gyj to characterize the strength of the 
gravitational wave field Q], Q], Q]. Qg^ is related to the power spectrum Ph{k) by the relation 

f^5»(^) = Y(^) Ph{k) (29) 

where kn = 2TiHo/c, and Ho is the present day Hubble constant. The quantity ilg^ is the present 
day ratio of energy density of gravitational waves (per unit logarithmic interval in k) to the critical 
density of the Universe pcrit — 

For simplicity, below we shall assume a simple power law behaviour for density parameter 
Qgw = ^gw{ko) ■ (k/ko)'^'^ . Although restricted, this form of spectrum is a good approximation for 
a large variety of models in gravitational wave frequency range of our interest. For example, this 
type of a power spectrum arises due to the evolution of relic gravitational waves with a (primordial) 
spectral index equal to ny, (i.e. Ph{k)\prim oc A;"^). The flat, scale invariant power spectrum (also 
known as Harrison-Zeldovich power spectrum) corresponds to ut = 0. In general the power law 
spectrum for Qgyj just assumes the absence of features in the spectrum of gravitational waves at 
the wavelengths of our interest. 

Let us consider electromagnetic radiation from a distant quasar. Expression (I24bl) allows us 
to calculate the expected angular fluctuation in the position of this quasar caused by a stochastic 
background of gravitational waves. In order to proceed we require to specify the limits of integration 
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kmin and kmax in (|24bp . kmin and kmax determine the frequency range of gravitational waves that 
can be probed by precision interferometry. The lower limit kmin is determined by the time duration 
of observations T^bs, kmin ~ ^"^/cTobs- The upper limit kmax ~ 2t: /c5t is determined by the time 
resolution of the observations 5t, and we shall assume 5t <^ Tots (i-e. kmax ^ kmin)- Let D be the 
distance to the quasar, which we shall assume is comparable to the Hubble length, i.e. D ~ Lh = 
cH~^. We shall be working under the assumption kD = InD/Xgw ^ 1, corresponding to the 
reasonable condition that the gravitational waves of our interest have a wavelengths much shorter 
than Lh- 

As can be seen from expression (\27h the behaviour of the transfer function Aa^(/c) depends on 
value of the quantity bire^kLn- In order to analyze the various possibilities let us introduce 



— {BTrkminLn) 



-1/3 



2.3 • 10" 



(30) 



where we have assumed kmin = ^n/cTobs, and T^bs = 10 yrs. In the above expression, and elsewhere 
below we set Ho = 75 ^/Mpc for numerical evaluations. 

For the angular displacement < Aa^ > , in the case e <C , substituting ([281 into ()24bp , taking 
into account the definition ()29p and integrating in the limits from kmin to kmax we get 



< Aa^ > 



gw \'^oj 



407r2 (l-nr/2) 



kn 



kr, 



3 ^gw (,kmin) ( Tobs 



407r2 (l-rar/2) \ Th 
In the opposite case of e » e* we get 



kn 

kmax J \ 

for e <C 



(l + z)2 + l 



(31) 



(1 - ut) 

2)^gyj {kmin) 



kn 



Tobs 



ko 



nx 



ku 
k 



kn 



[l + z- Vl + z) , 



ko 

for e ^> e 



[l + z- ViT^) 



(32) 



(1 - ut) V ^-ff 

In expressions ()3ip and (j32p Tfj = H^^ is the Hubble time and kn = 2-k /Lh- We have assumed 
kmax ^ kmin and used kn /kmin — Tobs/Tn- In the above expressions we restrict our analysis to the 
case riT < 1 which covers most of the practically interesting cases, including tlt = corresponding 
to a flat primordial spectrum of relic gravitational waves. In further evaluations below we shall 
set the redshift z = 4, which corresponds to a redshift with significant amount of quasar sources 
available for observations. 

The measurement of < Aq^ > for distant quasars by the planned interferometric projects [3], 



[13] would be able to constrain either Qgw or Qg^e^, depending on the value of e compared with 
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e^,. A null result in the measurement of < Aa^ >, in the case e <C e*, would place the following 
limit on the energy density of gravitational waves 0,gw (using expression (|31|) ) 

2' 



^gw < 4.0 • 10" 



(1 - nT/2) 



1 /ias 



10 yrs 



' obs 



(33) 



The above expression serves as the (weakest) upper limit on that can be set by precision 
interferometry measurements, irrespective of the value e (since, as will become clearer from the 
expression below, for values of e ^ e* this upper limit only becomes more stringent). In the case 



e ^ e*, from (j32]) . we get the following upper limit for the quantity ^.g^e^ 



^gwe^ < 3-7 • 10"^^ 



(1 - ut) 



1 /uas 



10 yrs 

Tobs 



(34) 



In expressions (j33p and (j34p we have introduced Aor-ms = V < Aa-^ >, which is the root mean 
square of the angular resolution of the interferometer. This precision is around 1 ^as for the cur- 
rently planned interferometers, and reaches 0.4 ^as for the proposed MAXIM x-ray interferometer. 
It is worth noting that, at this angular resolution, the upper limit on Q.gyj < 6.4 • 10~^ that can 
be achieved by precision interferometric measurements is comparable to the current limits set by 

LIGO 0, mv 

FigureSshow. .he con^aint, on energy density p^amete. and the yelocity patametet . 

achievable with an angular resolution of Aa^ms = 0.4 jias promised by the MAXIM project |17l |. 
The figure also shows the current constraints on the Op,,, parameter 271] , 2g] , 291] , 30(] , along with 

31C, y, y. 



sensitivity levels of some of the planned experiments 

An independent measurement of Vig^ (at a level below Vig^ < 10~^) by ground based inter- 
ferometers [3'], planned space borne interferometer LISA [26], or Cosmic Microwave Background 



anisotropy and polarization measurements 



111 ] would allow to place direct constraints on velocity 



parameter e. In this case, using (j32p . we can calculate the upper limit on e achievable by precision 
interferometric measurements 



e < 3.3 • 10" 



(1 - ut) 



10 



-10 



gw 



( AOr 



\ 1 /ias 



10 yrs 



f obs 



(35) 



This constraint corresponds to the region e > £=„ on Figure [3l 

In table U we summarize the predictions for r^^u,, and the corresponding upper limits on e, for 
some of the viable models that generate a considerable amount of stochastic gravitational wave 



backgrounds 



[24, 



12]. 
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FIG. 3: The upper limit on the energy density and velocity parameter e achievable by an interferometric 
observation of a source at redshift z — A with resolution Aa^ms = 0.4 fias and Tohs — 10 yrs. The shaded area 
shows the region that can be ruled out by precision interferometry together with other existing constraints. 
The horizontal lines on the graph show various constraints on the density parameter Qg^ (solid hairy lines 
- current upper bounds, dashed line - future sensitivity levels). 

V. CONCLUSION 

Although gravitational waves are yet to be detected, there are currently upper limits placed by 
observations. The CMB places the most stringent limits on gravitational waves of cosmological 
origin of Qgw < 10^^^ extrapolated to laboratory scale frequencies. It is crucial that this limit 
depends on extrapolation of data measured at the very long wavelengths (comparable to the Hub- 
ble length) down to wavelengths of the order of one light year and much lesser. Further more, 
these extrapolations assume a spectral index close to zero, i.e. a flat (scale invariant) spectrum of 
primordial gravitational waves. At frequencies relevant for precision interferometry observations, 
the most stringent constraint is placed by the Pulsar Timing measurements Qg^ < 2 x 10~^ at 



frequencies 10 ^ — 10 ^ Hz Q, and LIGO results Qgw < 6.5 x 10 ^ at frequencies 51-150 Hz 
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Theoretical Model 


Predicted f2g^ 


Upper Limit 


on e 


Upper Limit 


on 


nig 


Relic gravitational waves, ut 


= 






i 10-14 


e ^ 0.4 




mg < 5.2 


10 


-23 


eV 


Relic gravitational waves, ut 


= 0.2 aXu^ 


1 Hz 


^gw 


i 10-10 


e < 1.8 • 10 


-2 


TUg < 2.3 


10 


-24 


eV 


Local Strings 






^gw ^ 


« 10-8 


e < 3.9 • 10 


-3 


nig < 5.1 


10 


-25 


eV 


Global strings 






^gw ^ 


i 10-12 


e < 8.4 • 10 


-2 


mg < 1.1 


10 


-23 


eV 


Extended Inflation 






^gw ^ 


« 10-9 


e < 8.4 • 10 


-3 


mg < 1.1 


10 


-24 


eV 


1"* order EW transitions 






^gw ^ 


« 10-9 


e < 8.4 • 10 


-2 


mg < 1.1 


10 


-24 


eV 



TABLE L An upper limit on the velocity parameter e, for some viable models that predict a considerable 
gravitational wave background, that can be placed by interferometric observations with resolution Aarms = 
0.4 /ias, and observation time Tobs = 10 yrs. The last column shows the upper limits on the mass of the 
graviton rUg (see equation (|17p ) for a fiducial wavelength Xg^^ « 1 ly. 



28l ]. As can be seen from Table HI there are a host of viable theoretical models that predict gravita- 
tional wave backgrounds above the sensitivity levels of planned experiments like Advanced LIGO, 
LISA and SKA-PTA. If these experiments detect a gravitational wave background, the precision 
interferometry observations would be able to place strong constraints on e < 10"^. Let us finally 
note that, these interferometric measurements would be directly sensitive to gravitational waves 
with frequencies Vg^ ^ 2TT/Tohg ~ 2 • IQ-^ Hz, which is the frequency region that would be probed 
by LISA and SKA-PTA. 

Along side the candidates for a statistically isotropic gravitational wave background, we can 
also expect a significant stochastic background of gravitational waves from galactic white dwarf 
binaries {4] which are expected to be a dominant contribution to "noise" in LISA. Consideration of 
anisotropic gravitational wave backgrounds requires an approach slightly different from considera- 
tions of Section IIIIl but it is reasonable to assume that these sources could also place considerable 
upper limits on e. 
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APPENDIX A: EVALUATION OF THE TRANSFER FUNCTION 

Let us evaluate the integral in expression (j26p in the physically interesting case when e — > and 
kD ^ 1. In this case the expression for the transfer function can be separated into two distinctive 
contributions 



(Al) 



where Aq!^^(/c) is the non-resonance contribution 



Aa%ji{k) = — 



2\3 



dfi (1-^2) 



-1 



sm 



1 

+16 



+1 



dfi (l — Ai 



2\3 



l-e+Afi 



sm 



and Aa^(A;) is the resonance contribution 



1 



2\3 



sm 



(l-6-/.)2 

{e(^-^-/-)} 

(1-6-^)2 



(1 



(A2) 



(A3) 



l-e-A^l 

The quantity A/u occurring in the limits of integration in the above expressions is fixed by the 
condition for the resonance to occur. This condition corresponds to the region, around fi = 1 — e, 
where the sine function undergoes a few oscillations. Thus A/i = NXg^/D = 2'KN/kD, where N is 
the number of oscillations of the sine function, around the point /i = 1 — e, included in evaluation 
of the resonance. The value of N is limited by the condition A^u = lirN/kD <C e, implying 
N <^ ekD/2TT. Since in all our considerations we assume e ^ 1, and e^kD S> 1, the condition 
imposed on N is consistent with an additional condition ^ 1 that we shall assume. 

When evaluating ()A2p . in the case of e ^ 0, we can neglect the second integral in comparison 
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with the first. In evaluation fo the remaining integral we set e = 0. Thus, we get 

1 



AaW^) ^ j^jdfi (l-^)(l + ^f sin2(^^(l-^) 

-1 
1 

= ^1 dfi + fif (l - cos {kD (1 - /.)) 

-1 
1 



(A4) 



where, assuming kD ^ 1, we have explicitly separate out the rapid oscillatory part and neglect it. 

In order to evaluate (|A3p . in the case of e — > and kD ^ 1, it is helpful to notice that the 
factor (l — fi"^^^ in the right side of (|A3p is a slowly varying function over the range of integration. 
Taking this factor (evaluated at ^ = 1 — e) outside the integral we get the following approximation 
for the resonance part of the transfer function 



dfj. 



sm 



1 



+Ntt 



e'kD 



dx- 



sm^ X 



-Nit 



1 



-^e-'kD \l-0\j^ 



1 



1 c. 

-ire^kD. 
4 



(A5) 



Finally, the total transfer function, given by the sum of the non-resonance (lA4h and resonance 
parts (lASp . has the following form 



Aa\k) 



1 

20 



1 + bne-^kD 



(A6) 



APPENDIX B: THE SURFING EFFECT IN THE PRESENCE OF COSMOLOGICAL 
EVOLUTION 

In Section [TI] when evaluating expression ([9]) for simplicity and clarity we had neglected cosmo- 
logical evolution of the gravitational wave amplitude, i.e. we set a{rj) = 1. When we incorporate the 
cosmological expansion, it no longer becomes possible to evaluate the integral in expression Q in 
terms of elementary functions. Although this complicates the calculational aspects of the problem, 
much of the physical aspects remain the same as discussed in Sections HI] and IIIII Introducing 

-f = k{l-e-ke'), (Bl) 
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the expression ^ (in the case of a matter dominated universe governed by the scale factor ([2])) 
can be rewritten in the following form 

Mv,x') = ^e'e'^p,kWo I dx^, (B2) 

liVo-D) 



The expression ([15]) for the angular displacement due to a gravitational wave modifies correspond- 
ingly to 



'yvo 



dx—^. 



(B3) 



At this point it is convenient to consider the non-resonance and resonance contributions sepa- 
rately. Let us consider each in turn, beginning with the non-resonance contribution. For this case, 
I'Ho ^ 7(^0 — -D) ^ 1, and the integral in expression ()B3p can be evaluated asymptotically to give 



1 



AuNR ^ -7;ho e'e pik Kf -fk-q^ { ^ + ^ ( ^ 



1 



liVo-D) 
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— ho e e Pik kyl e '^''^° 
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- I e 
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i-yD 



Vo 



„2nLH(i-(i+^)-^/^) _(l + ^)- 



l — e — Le' 



(B4) 



where we have used 7]o = 2Lh and expression (jlOp to relate D and redshift z. The key difference 
between ()B4p and (jlSp is the presence of the cosmological redshift factor (1+z). This is the reflection 
of the fact that the non-resonance part of Aa probes the gravitational wave field at emission when 
the field was a factor (1 -|- z) stronger than at present. The non-resonance contribution to the 
transfer function (j25p can be evaluated in a fashion similar to the evaluation of Aa^^ in Appendix 
lAl The result is given by 



^a%j,{k) = ^ I cin \Aa{kY 



_ 1 
=0 ~ 20 



(l + z)2 + l 



(B5) 



Let us now turn to the resonance contribution in the transfer function. This contribution can be 
estimated by setting 7(7/0 — -D) ^ 7?7o <C 1. In this limit we can approximate the expression (|B3P 
in the following way 



AttR 



-ho e'e^pik krf -ykrjl - 
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(B6) 
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FIG. 4: The integrand figuring in the evaluation of the transfer function Aa|j.(fc) for the case of resonance, 
i.e. fi ~ 1 — e. The dashed fine shows the approximation used in (|B7|) . while the solid line shows the integrand 
in the exact evaluation. The value of redshift is set to z = 4. The overall (common to both the curves) 
normalization has been chosen arbitrarily. 

Thus, in comparison with expression (jlSp taken in the limit 7Z) 0, the expression ()B6P has an 



extra factor \/l + z. Using this approximation we can estimate the resonance part of the transfer 
function as follows 



1-e+Afi 



dfi (1 — /U 



2\3 



il + z)TTe^kD 



sm 



(1 
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^ne'kLHil + z)^2--^==j, (B7) 

where the resonance integral has been evaluated in the same manner as in Appendix [XI Figure H] 
shows the comparison of the approximate integrand used in ()B7p and the exact integrand calculated 
numerically. As can be seen from the figure, expression (|B7p gives a good approximation to the 
exact result. 

Combining expressions (IB4p and ()B7p we get the expression for the total transfer function 
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